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Spin-glass (SG) states of the 3-dimensional Ising Edwards-Anderson model under a static 
magnetic field h are examined by means of the standard Monte Carlo simulation on the field- 
shift aging protocol at temperature T. For each process with {T;t^,h), iw being the waiting 
time before the field is switched on, we extract the dynamical crossover time, tcr(r; t„, ft). 
We have found a nice scaling relation between the two characteristic length scales which are 
properly determined from ta and tw and then are normalized by the static field crossover length 
introduced in the SG droplet theory. This scaling behavior implies the instability of the SG 
phase in the equilibrium limit even under an infinitesimal h. In comparison with this numerical 
result the field effect on real spin glasses is also discussed. 
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In spite of extensive studies in more than a decade 
the nature of ordering in spin glasses have remained un- 
settled.^ A most notable problem is on the stability of 
the equilibrium spin-glass (SG) phase in a static field 
h. The mean-field theory predicts its stability up to a 
certain critical magnitude of h,^ while by the droplet 
theory^"^ the SG phase is unstable even in an infinites- 
imal h. In our opinion, a fundamental difficulty of the 
SG studies lies on the following fact: spin dynamics gov- 
erned by thermally activated processes is so slow that 
the time window of experiments on real spin glasses, not 
to mention that of numerical simulations, is insufficient 
for us to distinguish whether a SG-state we are observ- 
ing is in an asymptotic regime close to equilibrium or far 
from it. In such circumstances, a possible strategy is a lit- 
eral numerical experiment on nonequilibrium dynamics 
of model spin glasses without referring to their equilib- 
rium properties which have been extensively looked for 
by some ingenious, but artificial methods. Following 
this strategy, we have been extensively studying aging 
phenomena in the three-dimensional (3D) Ising Edwards- 
Anderson (EA) model^^^^ including the comparison of 
the results with those measured in real spin glasses. 

In the present work we address the problem of the 
SG-phase stability in a field by the strategy mentioned 
above, concentrating on the following field-shift (/i-shift) 
aging protocol. A system is instantaneously quenched 
at t = from T = oo to a temperature T below the 
transition temperature Tc, it is aged in zero field until 
t — when h is switched on, and the induced mag- 
netization, called zero-field-cooled (ZFC) magnetization 
and denoted as MzFc{i'), is measured as a function of 
t' = t ~ tv,. At i < isothermal and isobaric aging un- 
der h — proceeds and the SG short-range order, or the 
mean size of SG domains, Rrit), grows a.s,^'-^'^'^^ 

RT{t) = bTt^l'^^\ (1) 

with hT and z{T) being constants. The characteris- 
tic time scale after the /i-shift most frequently exam- 



ined is the one where the logarithmic time derivative of 
MzFc(i'): i-e-, S[t') = d{Mzvc[t')/h)/d\nt' exhibits a 
peak. This behavior of S{t') is interpreted to represent 
a crossover from the isobaric aging state under /i = to 
that under h > Q. The peak position of S{t') is called 
here the crossover time and denoted as tcr, though it has 
been frequently called the effective waiting time. It is a 
function of T, and h. We then expect that in the time 
range t' < tcr, which we call the transient regime of h- 
shift processes, SG subdomains in local equilibrium with 
respect to /i > grow within each SG domains which 
were in local equilibrium with respect to ft, = at t = 
Furthermore we intuitively introduce the growth law of 
mean size of the subdomains in the transient regime as 

Compared with the growth law of Eq.(l), the second 
term with ax is added in the exponent. 

The last step of our analysis is to relate the two charac- 
teristic length scales above-introduced, i?w = Rrit = tw) 
and i?cr = Rti{t' = tcr)-, to the field crossover length, Lh 
introduced in the droplet theory.^ The latter separates 
characteristic behavior of droplet excitations in equilib- 
rium by their size L, such that it is dominated by the 
Zeeman energy (- hL'^/'^) for L > Lh and by the SG 
free energy gap (^ TL^) for L < Lh. Here d is the spa- 
tial dimension, T the stiffness constant of domain walls, 
and the gap exponent. Explicitly Lh is written as 

Lh - lTh-\ (3) 

where 6 = (c?/2 — ^)^^ and It is a constant. Surpris- 
ingly it is found that, with a proper choice of the values 
of gt, the scaled plot Ra/Lh vs R^i/Lh of all the data 
obtained for different T, and h lie on a single curve. 
The consequence of this scaling is clear; the /i-shift aging 
process is nothing but a dynamical crossover from the 
SG state in ft = to the paramagnetic state in ft > 0, or, 
the SG phase is unstable under ft > in the equilibrium 
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limit. Furthermore, when the simulated results are sim- 
ply extended to the time window of the experiment on a 
real Ising spin glass Feo.sMng.sTiOa our picture of the 
dynamical crossover turns out to be consistent with the 
experimental data which were interpreted as an evidence 
of the AT phase transition predicted by the mean-field 
theory. 

We carry out standard heat-bath Monte Carlo (MC) 
simulations on the 3D Ising EA model with interactions 
obeying a Gaussian distribution with mean zero and vari- 
ance J which plays a role of units of energy and tempera- 
ture (with fee = 1). The strength of field h is represented 
by its associated Zeeman energy. The time is measured 
in unit of the MC steps per spin (mcs). In the present 
work we use the values of parameters br and z{T) in 
Eq.(l) obtained previously;^ {T;bT,z{T)) = (0.4; 0.82, 
14.8), (0.5; 0.80, 11.8), (0.6; 0.78, 9.8), (0.7; 0.78, 8.7) 
and (0.8; 0.76, 7.9). Note that ~ 0.95^^'^^ for this SG 
model. In the temperature range T < 0.8 and in the time 
range t 10^ where we carry out simulations, Rrit) is 
less than a few lattice distances.^ Therefore we examine 
systems with a linear size of 24, but an average over at 
least 6400 samples with different sets of interactions is 
taken for evaluating physical quantities in each /i-shift 
process discussed below. 

Let us begin our discussions with the typical results of 
MzFc{t') and Mpcit'), the field-cooled (PC) magnetiza- 
tion defined by MzFc{t') with = 0, at T = 0.6 shown 
in Fig. 1. Both exhibit a jump at t' = 1, which is natu- 
rally attributed to one-spin flips of individual spins whose 
internal field at t' = is smaller than and opposite to h. 
The number of such spins decreases in the ZFC process 
at t < and so Mfc(1) > Mzfc(1)- After the jump the 
system ages to an isobaric aging state under h gradually 
and irreversibly. Associated with this, subdomains in lo- 
cal equilibrium with respect to h are expected to grow 
obeying Eq. (2) , though we have not succeeded to observe 
them directly in simulation yet. Also shown in the figure 
is M{t') which is observed in h = 0.75 after shift-down 
from the FC process with h = 1.0 of a period t„. One 
can see that the magnetizations under h = 0.75 (1.0) 




hlOO 



1 10 100 1000 10000 

f 

Fig. 1. MpcW (f = t), AfzFc(t') with = 4096, and M(t') 
under h = 0.75 after the shift from h = 1.0 (tw = 4096). 
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Fig. 2. S{t') for the /i-shift processes at T = 0.6 with tw = 4096 
plotted vs Inf. The data are for h =0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 
0.6, 0.75, 1.0 and 2.0 from right to left. 

at t' ^ 10^ (10^) are independent of the history of the 
system at t' < 0, i.e., the system reaches the equilibrium 
paramagnetic state at this time range. For a small h such 
as /i = 0.1, on the other hand, we cannot see even a pre- 
cursor of saturation in Mpcit')- Our interest is therefore 
in behavior of MzFc(i') for 0.1 < h < 0.75. 

In Fig. 2 we demonstrate a set of S{t') observed in the 
/i-shift processes at T = 0.6 for various h but with a fixed 
Except for h = 2.0, S{t') exhibits a peak at t' = tcr- 
All tcr obtained for various sets of h and are shown 
in the inset of Fig. 3 below. For smaller fields, h < 0.1, 
we obtain tcr — tw, as has been observed since the very 
beginning of the SG aging study. For larger h, however, 
tci becomes significantly smaller than t^- 

As noted at the beginning of this letter, a key quan- 
tity in the droplet theory on the SG aging in a field is 
the field crossover length, Lh, and its role on the dynam- 
ical crossover of present interest becomes apparent if we 
compare it with the characteristic length scales, Ra and 
of each /i-shift process. As already defined, the latter 
are a function of T, and h. In processes with such a 
large h or/and that Ry, ^ Lh holds, the aging dynam- 
ics after the /i-shift is dominated by the Zeeman energy 
and so R^ — Lh is expected irrespectively of t^- In op- 
posite processes with Ry, -C Lh, we expect i?cr — Rw 
These qualitative features among R^,Rcy, and Lh are 
made more transparent by the scaling plots of Rcr/Lh vs 
R^/Lh- Even if we simply put i?tr(i') = RT{t'),^° i-e., 
ttT = in Eq.(2), we obtain the scaling plot, from which 
we can see the tendency of saturation of R^/ Lh at large 
Rvi/Lh- 

More interestingly, as shown in Fig. 3, all the sets of 
{Rcr/Lh, Rw/Lh) obtained in the /i-shift processes at dif- 
ferent temperatures, T = 0.4,0.5,0.6,0.7, and 0.8, with 
different sets of and h, lie on a single curve, if the val- 
ues of flT and It are properly adjusted as (T; ay, It) = 
(0.4; 0.06, 0.85), (0.5; 0.075, 0.83), (0.6; 0.09, 0.81), (0.7; 
0.115, 0.78) and (0.8; 0.125, 0.75). We emphasize here 
that the values of the parameters other than ay and 
It are those previously determined independently of the 
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Fig. 3. Scaling plots of i?cr / L/i vs Rw / i/j of the /i-shift processes 
at T = 0.4 ~ 0.8 with ay properly chosen. The line represents 
2/ = In the inset we plot tcr vs. tw in the /i-shift processes at 
T = 0.6 with h =0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, and 0.75 from 
top to bottom. 

present analysis; B = 0.20,^'^^ and so (5 ~ 0.77 in Eq. (3) 
with d = 3, and (6t, 1/2;(T)) already listed above. The 
adjustment of ot is essential to obtain a unique scaling 
curve for the data set at each temperature. The proper 
choice of It, which commonly scales both axes at each 
temperature, makes the scaling curves at different tem- 
peratures to lie top on each other. Physically, the main 
T-dependence of It can be attributed to that of the stiff- 
ness constant T which is expected to vanish at T^. 

The scaling behavior we have found implies that, for 
h of any strength, if a system is aged under ft, = up to 
for which Lh is satisfied, the aging dynamics 

after the /i-shift is governed dominantly by the Zeeman 
energy, i.e., the system exhibits paramagnetic behavior 
at t' ^ tcr- If iw is smaller than such a value, say = 0, 
it is faster for the system to reach the paramagnetic state 
as demonstrated in Fig. 1. 

To reach the nice scaling shown in Fig. 3, the term 
intuitively added in the exponent in Eq.(2) plays an es- 
sential role, though the scaling behavior is not so sensi- 
tive on the precise value of ar, i-e., even if it is changed 
by ±0.01 from the value listed above, the quality of the 
scaling does not become significantly worse. A positive 
Ot implies that the subdomains in the transient regime 
grow more slowly than the domains do in the isobaric 
aging, i.e., i?tr(^ = < Rrit)- We consider it quite 
plausible that the domain growth in h starting from a 
state already aged \n h' ^ h (including h' = 0) is slower 
than that starting from a state not aged at all. In this 
context it is worth noting that .Joh et al. explained their 
experimental result tcr < for relatively large h in terms 
of the reduction of the barrier energy for excitations due 
to the Zeeman energy. Numerically we have also ob- 
served that tcr < tw, but our interpretation for it is quite 
different from theirs. 

Our MC simulation alone cannot exclude a possible 
existence of the SG phase under sufficiently small fields 
{h < 0.1). However, the deviation of i?cr from R^, is 
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Fig. 4. x"{i^; t) aXT = 0.6 for = 2-k/ui = 64, tw = 8192, and 
h = 0.3 (filled circles). When its branch at t > tw is moved to 
the left by the amount tmv = tw — tcr = 5600 (open circles) 
it merges to the reference curve oi h = 0.3 (solid line) at t' = 
*mrg — 2000 ~ tcr (scc Rcf. 12 for further details of determining 
tmrg). The dashed line represents the reference curve oi h = 0. 
Those for h = 0.05, 0.1, 0.2 (not shown) almost coincide with the 
reference curve of /i = 0. 

clearly recognized even at small x = R^j/Lh- It originates 
from the fact that the excitation gaps of droplets with 
a fixed size distribute continuously and a finite weight 
exists even at zero energy. This point has been recently 
emphasized by Jonsson et aP^ who have called the con- 
sequence of these excitations the weak-chaos effect. Our 
data of Ra/Lh vs. R^/Lh at small, though not suffi- 
ciently small, X are consistent with their argument. We 
may better call this the weak-paramagnetic effect. 

One more comment is in order on the out-of phase 
component of the ac susceptibilities, x"('^;i); measured 
in the /i-shift aging protocol. It is evaluated via the 
fluctuation-dissipation theorem, and its typical behav- 
ior is shown in Fig. 4. In the isobaric regime at t < t^i 
X"(w;t) is rather insensitive to h as seen in the figure. 
Just after the /i-shift, however, it exhibits a sharp jump- 
up, and subsequently it tends to merge to the one ob- 
served in the isobaric aging under h from t = Q, which 
we call the reference curve. Qualitatively similar behav- 
ior has been observed experimentally.^^ According to the 
droplet theory, the aging part of x"{ijj; t) in the isobaric 
regime is given by a function of L^/Rrit), where La; 
is a mean size of droplets which can respond to the ac 
field of frequency u).^'^° A naive extension of this idea 
is that x"(w;t) after the /i-shift is dominantly governed 
by i?tr(i')' tfie mean size of subdomains which is smaller 
than i?^. In fact, as demonstrated in Fig. 4, the time 
scale that the x"(a;;t) merges to the reference curve well 
coincides with tcr, the peak position of S{t'). Thus the 
following picture on the aging dynamics around t' ~ tcr 
is extracted. At t' < tcr subdomains of mean size Rtr{t') 
grow up obeying Eq.(2), and at t' > tcr the subdomains 
become main domains in the system. For a process with 
Rcr — Lh the system is regarded to be already in a para- 
magnetic state and SG domains do not grow any more. 
For a process with R^ < L^, on the other hand, the 
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Fig. 5. Dynamical crossover field ha (sofid circles) due to Eq.(4) 
and the AT critical points (open circles) determined exper- 
imentally in [16] for Feo.sMno.sTiOa under fields. Shown in 
the inset are the same data plotted w.r.t. (Tc — T')/Tc double- 
logarithmically. The line represents the AT critical line. 



mean size of SG domains grows following Eq.(l) with a 
shift of the time origin by a certain amount. This growth 
is observed through physical quantities, such as x"{^]t) 
in Fig. 4, which associate a length scale smaller than 
the mean size of domains at t' . The value of MzFc(t') 
at t' ~ tci is, however, still significantly smaller than 
Mpc{t'), and it needs further long time for Mzfc(^') 
to catch up M-pc{t'), or the mean size of SG domains 
reaches Lh- 

So far described are the results of our simulation on the 
EA Ising SG model. In order to look for their relevance to 
phenomena in real Ising spin glasses, let us introduce the 
dynamical crossover condition defined by R„{f") / Lh = 
Xc with Xc being a constant nearly equal to unity. It is 
rewritten as 

hor = CTjtZ/to)-'/''^^"\ (4) 

with ct = {xclT/bTY^^ and to being the microscopic 
time unit. This equation gives us an estimate of the value 
hci, Tcr or when the values of the other two param- 
eters are specified, and when h ^ { <, )hci, ior example, 
the system is regarded as in a paramagnetic (SG) state. 

Katori and Ito (KI)^^ determined the critical temper- 
ature. Tat, of Ising spin glass Feo.sMno.sTiOa under h 
as the one where the Mzfc starts to deviate from Mpc 
within the observation time of iob — 10^ s at each tem- 
perature. Regarding this as the dynamical crossover with 
tw = *ob and Tat = Tcr (= T), we evaluate her by Eq.(4) 
with Tc = 21.5 K and to = 10~^^ s. This means that we 
simply extend the result extracted from our simulation 
to the time range of 10^** mcs as we did with success in 
our previous work.^^ Surprisingly, as shown in Fig. 5, the 
her — T plot thus evaluated with a;c — 1.5 coincides very 
well with the h — Tp^ plot that KI regard the AT crit- 
ical line.^ Similarly, the freezing temperatures T{{uj,h) 
extracted by Mattsson et al?^ through the ac suscepti- 
bility measurement are also satisfactorily reproduced by 
our analysis with Tcr = Tf, = 27r/a; and her = h. 
These agreements, whose details will be reported else- 
where, are quite promising for us to consider that the 
field effect also on real spin glasses is well described by 



our dynamical crossover picture. For h of the magni- 
tude less than 1 Oe, however, becomes astronomi- 
cally large, which explains, we consider, why MzFc/h 
measured in laboratory-time deviate from M-pc/h, and 
Mpc /h itself does so from the paramagnetic susceptibil- 
ity at T ~ Tc for small h.^^ 

To conclude we have demonstrated that the character- 
istic length scales associated with ft-shift aging protocols 
in the 3D EA Ising SG model exhibit unique scaling be- 
havior. Because of the latter, we have succeeded to reach 
one of the most important equilibrium properties of Ising 
spin glasses via simulations on nonequilibrium (aging) 
dynamics: their SG phase is unstable under a finite field 
in equilibrium. 
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Institute for Solid State Physics, the University of Tokyo. 
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